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Abstract. Kinetic Monte Carlo simulations are performed to study the influence 
of external magnetic fields on the growth of magnetic fee binary alloy nanoclusters 
with perpendicular magnetic anisotropy. The underlying kinetic model is designed to 
describe essential structural and magnetic properties of CoPt3-type clusters grown on a 
weakly interacting substrate through molecular beam epitaxy. The results suggest that 
perpendicular magnetic anisotropy can be enhanced when the field is applied during 
growth. For equilibrium bulk systems a significant shift of the onset temperature for 
LI2 ordering is found, in agreement with predictions from Landau theory. Stronger 
field induced effects can be expected for magnetic fcc-alloys undergoing L1q ordering. 
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1. Introduction 

It has recently been shown that CoPt3 ultrathin films and nanoclusters, grown by 
molecular beam epitaxy (MBE), display perpendicular magnetic anisotropy (PMA) 
with potential applications in high density magnetic storage media [U El • Previous 
simulations of CoPt3 nanocluster growth on a weakly interacting substrate suggested 
that PMA originates from the combined effects of Pt surface segregation and cluster 
shape |U E]. As suggested by experiments on Co-Pt multilayer systems jE|, out-of- 
plane CoPt bonds contribute to a magnetic anisotropy energy that tends to align the 
Co-moments along the bond direction and thus favours PMA. When a perpendicular 
magnetization is induced by a magnetic field applied during growth, this type of a local 
magnetic anisotropy energy can be expected to lead to an additional preference of out-of- 
plane relative to in-plane nearest neighbour Co-Pt pairs. MBE-growth of CoPt3-clusters 
in a perpendicular field should therefore improve the formation of PMA in nanoclusters. 

In order to demonstrate this latter effect we perform kinetic Monte Carlo 
simulations of growth within a binary alloy model, supplemented by a magnetic 
anisotropy energy H\ based on Co-Pt bond contributions. Isotropic magnetic 
interactions, on the other hand, do not contribute in a significant way to any growth- 
induced structural anisotropy. Their main effect in growth simulations is a small 
renormalization of chemical interactions, which can be neglected. However, both types of 
magnetic interactions in Co-Pt and related fcc-alloys can lead to an intriguing interplay 
between structural ordering phenomena and magnetism in samples at equilibrium [Jj. 
We confirm this by comparing predictions of Landau theory with equilibrium Monte 
Carlo simulations. 

2. Definition of the model 

The model we consider refers to binary fcc-alloys with composition AB3. A minimal 
set of chemical interactions between the atoms is chosen that is consistent with 
the main processes during structure formation: atomic migration in different local 
environments, surface segregation of the majority atoms and ordering with LI2- 
symmetry. On a semiquantitative level, these features can be reproduced by effective 
chemical interactions Vaa, Vab and Vbb acting between nearest neighbour A or 5-atoms 
on an fee lattice. The linear combinations I = {Vaa + Vbb — 2Vab)/4, h = Vbb — Vaa 
and Vo = (Vaa + Vbb)/2 determine the bulk ordering temperature T = 1.831 /ks [Hj, 
the degree of surface segregation of one atomic species and the average atomic binding 
energy, respectively. Following previous work |4 , we adapt our model to CoPt 3 (A = 
Co, B = Pt), where T ~ 960 K, and Pt-surface segregation is strong (nearly 100 % 
jU]). Compatible parameters are / = 1, h ~ 4 and Vq ~ —5, where we have used 
fcsTo/1.83 ~ 45meV, corresponding to 523 K, as our energy unit. The value for Vo 
describes the average binding energy for intermediate coordination numbers experienced 
by atoms within the growth zone near the cluster surface. The substrate (111) surface 
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is weakly attractive with a potential V s = —5 that acts on adatoms in the first layer. 
Compared to a Pt(lll) surface with three bonds of typical strength Vo, this amounts 
to about 1/3 of the energy of a single Pt-Pt bond. 

The elementary processes in our continuous time kinetic Monte Carlo algorithm, 
which drive the cluster growth, are (i) codeposition of A- and I?-atoms with ratio 1:3 
and total flux F, (ii) hopping of atoms to vacant nearest neighbour sites, and (Hi) 
direct exchange of unlike nearest neighbour atoms; one of them is an adatom with low 
coordination (3 to 5) on top of a terrace and the other one a highly coordinated atom 
(coordination 8 to 10) underneath. Such exchange processes facilitate Pt segregation to 
the surface, which otherwise becomes kinetically hindered through the incoming flux. 
In our simulations we choose F = 3.5 monolayers per second. As far as possible, other 
kinetic parameters are adapted from known diffusion data: Jump rates for A and B 
atoms are of the form uexp[—(U + max(0, &E))/k&T] where v ~ 8.3 x 10 11 s _1 is the 
attempt frequency, U ~ 5 jTU] is the diffusion barrier, and AE denotes the energy 
difference before and after the jump. Direct exchange processes among pairs of unlike 
atoms are subjected to an increased barrier U + U x with U x = 5. For more details on 
that model and on the choice of these parameters we refer to [3]. 

As explained in the introduction, PMA in CoPt3-nanoclusters can be related to a 
structural anisotropy, expressed by the parameter 

P = (n ( l oPt -n^ oPt )/N, (1) 

which is defined in terms of the difference between numbers of Co-Pt bonds out of plane, 
n^ oPt , and in plane, n^ 01 ^. N is the total number of atoms in the cluster. The underlying 
analysis [3] can in particular account for the temperature window where PMA occurs, 
which lies below the onset of Ll 2 ordering. Furthermore, it provides an interpretation 
of PMA in terms of cluster size and shape, and the degree of Pt-surface segregation. 
The associated magnetic model assumes a local crystalline anisotropy energy that 
involves only Co-moments fi as its dominant part. Within a bond picture, each Co- atom 
experiences an anisotropy energy that is a sum over Co-Pt bond contributions 

-A(».6r/(\»\\d\y. (2) 

Here d is a bond vector connecting the Co-atom with one of its nearest neighbour 
Pt atoms, and A ~ 225 /xeV is a parameter deduced from the interfacial part of 
the measured magnetic anisotropy of Co-Pt multilayers jHJ- We remark that these 
measurements also appear to justify the bond model implied by (j2J), when two different 
interfacial orientations are compared. For (111) interfaces the anisotropy energy per 
surface area is found about twice as large as for the (100) orientation. Considering 
the different angles of bonds to the surface and the different packings, equation (J2J) 
indeed can reproduce this difference with one consistent value for A. Since A > 0, 
spin alignment along a Co-Pt bond is favoured. Certainly, A is small relative to the 
chemical interactions Vaa-, Vab and Vbb as we U as I an d h. Hence, in the absence of a 
magnetic field it is justified in a first approximation to neglect (J2J) in simulating growth, 
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but to include it a posteriori in order to relate a given cluster structure to its magnetic 
properties. 

3. Growth in external fields 

However, when an external magnetic field is applied, we show here that inclusion of 
the magnetic anisotropy energy (J2J) in the growth process itself leads to a small but 
notable change in the clusters' short range order such that PMA gets enhanced To 
demonstrate this, we apply a strong field B s perpendicular to the substrate that drives 
the magnetization towards saturation. Co-moments /j, in (J2J) are then aligned along the 
[111] direction, which introduces an asymmetry in the probabilities for jumps that form 
or break CoPt-bonds. For in-plane bonds, (J2J) gives no contribution, whereas the energy 
of an out-of-plane Co-Pt bond is changed from Vab to V' AB = Vab — \ A in the fee lattice 
geometry. 

Clearly, the expected magnetic field effects on growth are small in view of the 
smallness of the magnetic anisotropy parameter A in (J2J. To obtain reliable results for 
the change in the structural anisotropy parameter P with and without the field, we have 
carried out simulations including several much larger parameters A up to A = 1. Note 
that the energy unit chosen (see above) implies that the physical value A ~ 250/ieV 
corresponds to A = 5 x 1CT 3 . Results for those different artificial A-values will be used 
in turn to extract the physical change of P by interpolation. 

From this analysis we first show results for the dependence P = P(N, A) on cluster 
size N, see figure 1. All data sets referring to different A exhibit an approximate linearity 
in iV" 1 / 3 , 

P(N,A) ~Ps(A) N~ 1/3 . (3) 

Concerning the iV-dependence, this confirms that the structural anisotropy essentially is 
a surface effect!" . The figure includes the case A = 0, represented by open symbols. This 
data formally corresponds to zero magnetic field, because for B = at the temperatures 
considered the Co-moments \x are oriented randomly so that (J2J) gives no contribution 
to P. The important observation is that P(N, A) for fixed N significantly increases with 
A. Figure 2 explicitly shows that the behaviour of the coefficient p s {A) obtained from 
fits according to Q, appears to be consistent with a linear A-dependence. Assuming 
linearity even in the limit A — > 0, we obtain (p s {A) — p s (0))/p s (0) ~ 3 • 10~ 2 for the 
physical value A = 5 x 10~ 3 . 

This result immediately translates to the magnetic field-dependence of the 
structural part of the total magnetic anisotropy energy E s , which is proportional to 
NP. Using (PI), we can write 

E s = K S (B)N 2 / 3 , (4) 

J Fits of data in figure 1 by straight lines and extrapolation to TV -1 / 3 — > yield in addition a non-zero 
bulk contribution to P which, however, is negligibly small in comparison with P-values for clusters 
with N ~ 10 3 atoms. 
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Figure 1. Simulated structural anisotropy parameter P depending on cluster size 
N in the presence of the saturation field B s at temperature fcsT = 1.2 (T ~ 630K). 
Averages are performed over 20 clusters. Full lines display linearity in TV -1 / 3 . Open 
symbols are equivalent to the case of zero magnetic field. 



< 

D. 



6 
5.5 

5 
4.5 

4 
3.5 

3 
2.5 

2 

1.5 4 
1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
A 



Figure 2. Coefficient p s {A) determining the structural anisotropy versus the magnetic 
anisotropy constant A at T = 1.2. Data points result from the slopes of figure 1 
including the case A = 0.05. Dotted line: linear fit. 

again with the relative change (K S (B S ) — K s (0))/K s (0) ~ 3 • 1CT 2 . Thus, our main 
result is that the surface anisotropy constant K s in (HJ) increases by about 3 % when a 
perpendicular field of the strength of the saturation field B s is switched on. Near T ~ 1 
(T ~ 523i^), where the anisotropy develops its maximum JT], we obtain analogous 
results (not shown). 
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4. Interplay of structural and magnetic order in the bulk 

So far we have shown that in the presence of a strong perpendicular magnetic field the 
anisotropy energy if a notably affects the frozen-in structure of growing nanoclusters and 
thus improves PMA. Based on Landau theory, we now study effects of external magnetic 
fields on the equilibrium structure of bulk systems or homogeneous films. Effects of this 
type are known in principle but have been explored only for few specific materials [T2"j- 
[T3*] . For CoPt3, quantitative calculations have been performed in the past within the 
cluster variation method in the tetrahedron approximation Landau theory, which 
we use here allows us to incorporate the magnetic anisotropy in a simple manner and 
helps clarifying the interrelation between different effects. 

For the free energy / per Co-atom in CoPt3 in the presence of an external field B 
we propose the form 

f = fs + f M + /int - M • B, (5) 

where M = (Mi, M2, M3) is the average magnetic moment per CoPt3-unit. The first 
term in (jSJ denotes the structural free energy that describes Li2-ordering without 
magnetic contributions [T5] . 

3 / 3 \ ^ 

/^Ef^+i^ + i +™hM*- (6) 

a=l ^ ' \a=l / 

The structural order parameter components ipa, at = 1,2,3, describe amplitudes for 
layering along one of the cubic axes a of the fcc-lattice such that CoPt-layers and Pt- 
layers alternate. Superposition of these layer structures in all three directions yields the 
Ll 2 -structure. As usual, we set r(T) ~ r (T — T sp ); r > 0, u > 0, < v < u and w > 0, 
where T sp denotes the spinodal temperature. The associated ordering temperature To 
is determined by r(T ) = (2/9)u> 2 /(3w + v). 

The second term in (jSJ) is the isotropic magnetic free energy, with the behaviour 

/ M (M) ~ ^M 2 , (7) 

as M — > 0, where b(T) ~ b (T — T c ) and T c denotes the Curie-temperature of a 
disordered alloy. Finally the third term in (j3J), / int describes the coupling between 
structural and magnetic order. Retaining only the lowest-order terms allowed by 
symmetry, we write 

/int = |fE^) M2 + |E^ ( 8 ) 
\a=l J a=l 

For later purposes we define the Landau coefficients such that ip a becomes dimensionless, 
with \ip a \ = V'max = 1 for perfect structural order. Expressions (0) - © imply far- 
reaching interdependencies between structural and magnetic properties. Combination 
with simple mean-field arguments allows us to make order-of-magnitude predictions 
which can be tested experimentally. 
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First, collecting terms in (JSJ which are quadratic in ip a , we find that in a 
state magnetized parallel to the [111] direction, both the spinodal and the transition 
temperature for Ll 2 -ordering are shifted by an amount T sp (M) — T sp ~ Tq(M) — Tq = 
AT (M) with 

ATq(M) = -(3d + c 2 )M 2 /3r . (9) 

To determine the c-coefficients in (JHJ) we note that the same combination 3ci + c 2 enters 
the difference in Curie temperatures between disordered and ordered samples, 

ATcOmax) = T C (?/W) ~ T C = -(3ci + C 2 )/& - (10) 

This relation is immediately obtained from (JHJ) by requiring the terms quadratic in M 
to vanish. For the ordered samples we have set \if) a \ ~ VVax, because the discontinuity 
■00 in ipa at To is large and Curie temperatures are considerably lower than T ~ 960 K. 
Within the same consideration the inverse magnetic susceptibility x~ l develops a 
discontinuity at the order transition at T of magnitude (3ci + c 2 )-?/>o- Experimentally, 
Tb(VWx) — 400 K, while ATc(ip ma , x ) > 150 K [T3j. The reason for an increased Curie 
temperature in quenched, disordered samples is the increased number of Co-Co direct 
exchange couplings relative to the ordered state. This temperature shift already allows 
a rough estimate of 0- Suppose, for example, the magnetic field is strong enough to 
drive the magnetic moment M to values near its saturation value M s . Using mean 
field arguments, we expect that both r /k B and b M^/k B are of order unity. Therefore 
the two expressions |T0|) and (JHJ) should be of similar magnitude. In other words, for 
very strong magnetic fields B || [111] equation (JHJ) predicts a lowering of To by about 
AT 0imax ~10 2 Kin CoPt 3 . 

A further effect is worth mentioning. When the orientation of the external field 
is changed from [111] to [001] and c 2 < 0, the disordered phase develops an instability 
towards formation of a layered structure with ipi — ip2 — and ipz 0, i.e. with 
L1 Q symmetry. The corresponding transition point is increased relative to the spinodal 
temperature for [111] field orientation, T sp (M), by an amount (2/3)|c 2 |, but will exceed 
To(M) only if |c 2 | is sufficiently large. Otherwise this layered structure cannot form 
because the conventional ordering transition at T (M) will set in before. In CoPt3 
anisotropic magnetic interactions determining c 2 are weak (see below) so that this effect 
is unlikely to be observable. 

At this point let us attempt to explicitly relate the parameters c\ and c 2 to a specific 
magnetic model. For that purpose we adopt a phenomenological spin-Hamiltonian 
H M = i? ex + i?A pertaining to CoPt 3 alloys. As its main contribution we assume 
classical Heisenberg-type isotropic exchange interactions if ex within nearest-neighbour 
Co-Co and Co-Pt pairs. This simplified form seems sufficient as we are dealing only with 
the energetics of homogeneously magnetized states. The respective magnetic energies 
for Co-Co and Co-Pt pairs with parallel spins are — J Ml and — J'M 2 , respectively^ In 

§ Pt-moments have been shown to depend on the actual atomic environment |14j but this effect is 
neglected here. 
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Figure 3. Simulated L±2 order parameter W(T) for a bulk system at equilibrium as 
a function of temperature. Averages were performed over 6 realizations. Among the 
magnetic interactions only the anisotropy Ha has been included in the simulations. 
The upward shift of the ordering temperature Tq(B) — To ~ 2 K is clearly recognizable. 

addition, an anisotropy Ha connected with Co-Pt bonds, see (J2J), is included. Within 
first order perturbation theory the change in free energy due to Hm is given by an average 
of Hm over an unperturbed ensemble specified by the ip a . Mean-field calculations then 
yield 

ci = J - 2J', c 2 = -2A/M 2 S . (11) 

In order to test the validity of the Landau approach based on the coupling (JHJ) 
between structural and magnetic order and its prediction (0), we first focus on the 
anisotropy energy Ha (ignoring H ex ) and perform Monte Carlo simulations of the lattice 
model described before, but now under equilibrium conditions. The simulation box 
contains 43200 atoms and is subjected to periodic boundary conditions in all three 
directions. Equilibrated configurations are generated by allowing any pair of unlike 
atoms to exchange positions. In figure 3 the structural order parameter is plotted 
versus temperature, and indeed shows an increased transition temperature in the fully 
magnetized state in comparison to the zero field case. We have again used A ~ 225 /xeV, 
see above, which gives an increase of To by 2 K. This increase favourably agrees with 
(JOJ) when we set c\ = 0, see (jlljl . and r^/k-Q ~ 1. These results show that with realistic 
parameters for the microscopic anisotropy energy a measurable shift in the ordering 
temperature is observed. The physical reason for the increase of To in figure 3 is that 
through the presence of Ha the six nearest neighbour sites of a Co-atom with bond 
vectors not orthogonal to [111] become preferentially occupied by Pt-atoms during the 
equilibration process, thus favouring the ordered structure. Basically this is the same 
mechanism that yields an increased structural anisotropy parameter P discussed before. 
Figure 3, therefore, provides additional evidence, that Ha affects the ordering transition. 
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When exchange interactions are included, equation Q predicts a much larger 
downward shift in T , which we can estimate now in a more quantitative manner. From 
[Hj and known Co- and Pt-moments we obtain JM 2 ~ 140 K and J'M 2 ~ 25 K, while A 
is much smaller. Thus (3ci + C2)M 2 /3 ~ 90 K, in agreement with the order of magnitude 
estimate of (jHJ) given before, which was based on the comparison with measured Curie 
temperatures. 

Again, the corresponding change in T according to Q can be compared with 
predictions from numerical computation. However, within a lattice-model with nearest- 
neighbour chemical and isotropic exchange interactions, extra simulations for this 
problem actually are not needed when we note that in a fully magnetized sample the 
effect of H ex is simply equivalent to modified chemical interactions 

V AA = V AA -JM 2 ; V BB = V BB ; V' AB = V AB -J'M 2 S . (12) 

Since the ordering transition temperature satisfies /cb^o — 1-83(Vaa+Vbb — 2Fa_b)/4, we 
find, using (O, k B AT {M s ) ~ — 1.83( J — 2J')M 2 /A. Because of (dlj) and |c 2 | < c x this 
is consistent with (JHJ) and establishes a semi- quantitative estimate ATo(M s ) ~ — 41 K. 
Clearly, since the exchange terms JM 2 and J'M 2 in ()12j) are relatively small, we could 
neglect them in our previous discussion of structural anisotropy. 

Experimentally, for CoPt3 it is probably very difficult to induce a sufficiently 
strong magnetization near T such that the shift AT oc M 2 becomes measurable. The 
situation is more favourable for Llo-ordering of Co^Pti-^ alloys near x = 0.5, where 
the Curie temperature is closer to the ordering temperature T and even intersects 
with To on the Co-rich side in the phase diagram. A similar situation occurs for 
Ni-rich Ni^Pt^-i alloys [Zj. For Llo-ordering the structural part of the free energy 
fs is again based on ©, however, with a change in sign in w at x = 1/2, v 
sufficiently negative and additional stabilizing higher order terms ^B]. Elastic energy 
contributions due to the tetragonal distortion in the Ll phase are thereby neglected. 
Coupling between structural order parameters and the magnetization can again be 
represented by equation ((HJ). Considering magnetic fields B || [001] we now obtain 
AT sp (M) ~ AT (M) = -( Cl + c 2 )M 2 /r and AT c (Vw) = ~{c x + c 2 )/b instead of 
(jHJ) and with a measured value of about 10 2 K for the latter temperature shift. 

Using similar arguments as before, both temperature shifts AT sp (M s ) and AT c (^ max ) 
are expected to be of the same order. 

5. Conclusions and outlook 

Effects of external magnetic fields on both the far-from-equilibrium MBE growth of 
magnetic fee alloys and their bulk structural phase behaviour at equilibrium have been 
studied. The statistical model we employed pertains to binary alloys of the CoPt3 type 
and Ll 2 -ordering. KMC simulations are supplemented by equilibrium considerations 
based on Landau theory. Using realistic parameters for CoPt3 our main aim was 
to explore the influence of the local crystalline magnetic anisotropy energy, modelled 
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by a bond Hamiltonian Ha, on the structural properties of alloy nanoclusters, when 
a strong magnetic field is applied during growth. Such an effect, although small in 
CoPt3, appears to be general and indeed leads to a notable enhancement of the clusters 
structural anisotropy and the associated PMA. Experimentally, in order to induce a 
large magnetization, the substrate temperature should fall below the Curie temperature 
of disordered samples but at the same time remain in the known temperature window 
where PMA occurs. 

More favorable conditions for studying these effects experimentally should exist for 
alloys like CoPt or FePt [T7j undergoing Ll -ordering, if clusters could be grown along 
the c-axis. As known from thin film measurements [IB], the magnetic anisotropy in 
those alloys is a bulk property connected with alternating Co(Fe)- and Pt-rich layers 
caused by the Llo structural order. Hence, under growth conditions where Llo-ordering 
is suppressed, one can expect a substantially larger field-induced PMA. 
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